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A regularization for effective field theory with two propagating heavy particles is 
^ \ constructed. This regularization preserves the low-energy analytic structure, im- 
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1 Introduction 



One of the necessary ingredients in an effective field theory (EFT) is the 
ultraviolet regularization. A natural regularization for the EFT is one that 
removes from loop integrals all contributions with momenta above the EFT 
energy scale, Aeft, and that also implements a low-energy counting by making 
the small energy scales in the EFT, namely small momenta and light particle 
masses, the scales that set the magnitude of power ultraviolet divergences. One 
example of such a natural regularization is dimensional regularization in EFT 
with only light particles, such as Chiral Perturbation Theory (ChPT) [1,2] 
for quasi- Goldstone pions. In dimensional regularization the scale of ChPT is 
given in terms of the pion decay constant F^r, and power divergences of loop 
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integrals are proportional to the small scales (M^ and/or invariants formed 
with the small external momenta) divided by the appropriate powers of F^^. 
Although a momentum cutoff can be used, it does not deal well with power 
divergences, as these contain pieces proportional to powers of the cutoff. In 
this case, an EFT can also be formulated, but it is more comphcated to deal 
with because the low-energy power counting must be recovered in the process 
of renormalization. 

In an EFT with a heavy stable particle, such as ChPT with nucleons, the heavy 
particle mass M is a new scale in the theory (some heavy masses could be much 
larger than Aeft, e.g., baryon masses in large- A^c QCD). With the inclusion of 
M, dimensional regularization ceases to be a natural regularization. Indeed, 
the power divergences lead now to contributions proportional to fractional 
powers of M, i.e. contributions that give rise to terms that are powers of 
times logM^ in d — > 4 dimensions. One such example is the nucleon 
self-energy in ChPT calculated in dimensional regularization [3]. Although 
the EFT can be consistently regularized with dimensional regularization, a 
natural regularization is more convenient. 

In the case, where a single heavy particle propagates through a Feynman 
diagram, one way to remove fractional powers of M [7] is to perform an ex- 
pansion in 1/M at the level of the Lagrangian [4], eliminating in this way 
M from the propagators. This approach implies that all the Green functions 
must be analytic functions of 1/M in the low-energy region. This, however, 
is not correct because there are low-energy singularities in 1/M, such as the 
anomalous thresholds (e.g., in the scalar form factor of the nucleon [5]), which 
the EFT must retain. Recently, a regularization based on dimensional regu- 
larization that avoids these problems of the 1/Af expansion was developed. 
Building on work by Ellis and Tang [6], Becher and Leutwyler [7] developed a 
regularization scheme that is natural, preserves the low-energy structure and 
internal symmetries, and is Lorentz covariant. This regularization has been 
called infrared regularization. 

In the case of two heavy particles, the 1/M expansion cannot be applied due 
to the two-heavy particle threshold in the s-channel. It is, therefore, very im- 
portant to develop a natural regularization for this case as well. In this letter 
such a regularization is developed along lines similar to those of the infrared 
regularization [7]. Starting from standard dimensional regularization, a mod- 
ification of the Feynman parameter integration domain is introduced leading 
to a natural regularization. Its implementation is presented in a Lorentz co- 
variant framework for basic one-loop diagrams. 
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2 Reguleirization of one-loop diagrams 



In this section, the basic one-loop diagrams in the two-heavy-particle sector are 
discussed: these are the bubble diagram Fig. 1, the triangle diagram Fig. 2 and 
the box diagram Fig. 4. In each case, the dimensionally regularized Feynman 
integral / can be decomposed into 

I^I + R, (1) 

where the fractional powers of M are entirely contained in R. In addition, R 
is an analytic function in the low-energy domain (i.e. local) and is referred to 
as the regular part. The low-energy part / contains all the low-energy analytic 
structure of the integral. Because / does not contain any fractional powers of 
M, it satisfies a low-energy power counting as discussed in Sect. 3. The loop 
integral in the EFT Dimensional Regularization (EFTDR) is then identified 
with /. 

The low-energy domain is expressed in terms of the heavy mass M, the light 
mass m, and the invariants Pj • Pj formed with the external heavy particle 
momenta, and is defined by the conditions | Pi ■ Pj — | <^ A|prp and m <^ 
Aeft- As in ChPT, low energy quantities such as m, external momenta of 
light particles, differences of external momenta of heavy particles, etc., will be 
said to be of order p. 

Although integrals that carry factors of the loop momentum in the numerator 
are not explicitly presented here, such integrals can be reduced to combinations 
of scalar integrals of the form discussed here, to integrals that involve only one 
heavy propagator that can be regularized using the infrared regularization, and 
to integrals that involve only light propagators. 

2.1 The bubble diagram 




Fig. 1. 

The simplest diagram in the two-heavy-particle sector is the bubble diagram 
of Fig. 1 with the two heavy particles propagating in the s-channel. The scalar 
integral for this diagram is 
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where s = and P is the total incoming momentum. Combining the prop- 
agators using a Feynman parameter integral, the integration over momentum 
yields 



s 



d/2-2 



l2 = K2 dw -K - 1) + - ie , (3) 

-1 L4 J 



with 



This integral has a pinch singularity at w = when s — s+ = 4M^, which 
is the branch point of the two-heavy-particle cut. The integration boundaries 
give rise to fractional powers of M which spoil low-energy power counting. 
This can be shown below threshold by scaling the small quantity (M^ — s/4) 
out of Eq. (3): 



h = «2 
where 



(iY"'\''-^ I'dw (w^ + if-^ , (6) 



(6) 

is a small quantity in the low-energy domain which characterizes the infrared 
behavior of the bubble diagram. From Eq. (5) it is clear that only the integral 
limits can change the power of a. 

The integral can then be separated into its regular and low-energy part: 

72=l2 + i?2, (7) 

with 



oo 



7 2 = «2 (^) ' a"-^ j ^"dw {W^ + l)'^/2-2 _ (g) 

I2 and I2 have the same low-energy analytic structure as can be verified for 
d — > 4, where 
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167r2 



— 2 + 2q; arctan 



(9) 



and where, after performing one integration by parts, 



= (10) 

I2 is, therefore, given by the term of order a in the Taylor expansion of Eq. (9) 
in powers of ct, as can be immediately seen from 



q; arctan f—^ = —a — a arctan a, (11) 
\aj 2 

The term a arctan a is analytic in — s in the low-energy region. It is then 
clear that 1 2 and I2 have the same low energy analytic structure. 

In arbitrary d-dimensions, 1 2 is proportional to (s_|_ — s)('^~^)/^/y^, which 
imphes that it does not contain fractional powers of M. On the other hand, 
for d — > 4, i?2 is given by 



M''-^ i , , 

it2 = 2— — -K2 — ——^Oioxciajia^ (12) 
a — 3 Stt^ 

while in the low-energy domain and for arbitrary d 

i?2 = 2«2M^- E \4 . (13) 

where 



C = ^ ^2 tl ('14) 

((i-3-2n)m!n!r(f -l-m)r(f-l-n)' ^ ^ 

which shows that R2 is given by terms analytic in s — multiplied by a 
fractional power of M . In particular, this shows that the analytic terms not 
included in I2 are all of fractional-power type. 

The separation into a low-energy part and a regular part with the desired 
properties has been achieved. Indeed, from Eqs. (6) and (10) it is seen that 

I2 has two cuts: a cut due to the factor — s starting at s = s+ and a 
cut along the negative axis ending at s = due to the factor of 1/^/s. The 
discontinuity of I2 and I2 across the cut on the positive s-axis are equal, and 
therefore 1 2 reproduces the low energy analytic structure of I2. The cut along 
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the negative s-axis is located far outside the low-energy region. Moreover, I2 
does not carry any fractional power of M. On the other hand, R2 is purely of 
fractional power type, and has only the cut along the negative s-axis. 



2.2 The triangle diagram 




Pi 



Fig. 2. Triangle diagram in the two heavy-particle-sector. The solid line represents 
the heavy particle, the dashed line the light particle. 

The next diagram to be considered is the triangle diagram shown in Fig. 2 
with two heavy particles propagating in the s- channel and interacting via a 
contact interaction and the exchange of a light particle of mass m. 

The scalar integral for this diagram is 



Combining the heavy propagators with a symmetric Feynman parameter in- 



and subsequently combining this result with the light propagator using 




(15) 



tegral. 




(16) 




(17) 
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the momentum integration yields 



r / \ ](i/2-3 

h = 1^3 dz dw z C[w, z) — ie 
C(w, z) = {l- z)m^ + zM'^ - z\l - w^)^ 



z{l - z)-{Pl + P^') - wz{l - z)-{Pl - Pi 



In order to decompose 1^ into its low-energy and regular parts, 

J3=73 + i?3, (19) 

the integration range in w is first extended to the entire real axis in analogy 
to the bubble diagram, 



- rl roo r / \ l.d/2-3 

I3A = 1^3 dz dw z C(w, z) — ie 

JO J-oo L J 



(20) 



/3A contains the proper singularity structure associated with the two-heavy- 
particle threshold at s = s+ but does not contain the two-particle cuts cor- 



responding to the heavy-light thresholds at Pf = (M + my 



1,2). All 



this becomes clear by solving the corresponding Landau equation^ that de- 
termine the location of the singularities in the (w, 2)-parameter plane shown 
in Fig. 3. The two-particle thresholds are endpoint singularities, i.e., they lie 
on the boundary of the original integration domain [—1, 1] x [0, 1]. The two- 
heavy-particle cut (indicated by x in Fig. 3) occurs at (0, 1) and is unaffected 
by extending the domain as in Eq. (20). On the other hand, the heavy-light 



thresholds (indicated by •) are located at (+1, 



and 



-), and are 



lost. The heavy-light thresholds, indeed the full low-energy analytic structure, 
are restored by defining the low-energy part as 



-^3 — -^3A + -^3B , 



with 



- f°° f f°° \ I / -|rf/2-3 

/3B = — /«3 J dz ij dw + J dw j z C{w,z) — ie 



(21) 



f22) 



It is evident from Fig. 3 that I^, properly incorporates all two-particle thresh- 
olds at low energies and does not pick up any new low-energy endpoint singu- 
larities on the physical sheet .0 It is also important to note that J3 receives no 



^ For a review on Landau equations see Refs. [8-10] and references therein. 

In particular, the pseudo-thresholds = {M—my located at (±1, —m/{M—m)) 
and indicated by o in Fig. 3 are not promoted to thresholds on the physical sheet. 
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additional contribution from the leading (i.e. pinch) singularity (dark shaded 
area in Fig. 3) other than that already present in I^. This ensures that the 
analytic properties of and /a coincide in the low-energy domain and implies 
that is analytic there. 

By virtue of Eqs. (19) and (21) the regular part is 

i?3 = j'^d^z (^j dw + ^°dw^ z [C{w, z) - ie] ''^^"^ . (23) 

The dependence on fractional powers of M stems from the corner points lo- 
cated at (±1, 1) (indicated by ® in Fig. 3) and is entirely contained in i?3. It is 
instructive to verify this dependence on fractional powers of M explicitly. The 
denominator in Eq. (23) is non- vanishing and of zeroth order in the low energy 
expansion in the whole integration domain and can, therefore, be expanded 
in rn? and the small variables s+ — s, Pf — M^: 

C{w, z) = M'^z'^w'' + 0(s+ -s-Pf- M^; m') . (24) 

The leading contribution in an expansion of R3 is then 

fi3 = 2.,M"-^-^(l + ...) (25) 
11^2 _2_,_logi£l)(l + ...), (26) 



167r2 2M2 \4:-d ' 47r/i2 

where the omitted terms are analytic in the small quantities. The reason why 

i?3 has an ultraviolet divergence despite I3 being ultraviolet finite is that in 
order to contain the fractional powers of M while being regular in the low- 
energy domain, it must include the appropriate logarithmic dependence on /j,, 
and this can only emerge if there is a l/(d — 4) singularity. 

The low-energy parts are shown below for d — > 4. The integral I^a is given 
by: 



!3A 



1 /•! 1 

-r / ' (27) 

T\/S Jo y/Co 



with 

Co = zM^ + (1 - z)m^ - ^ - 1.(1 - z){P^ + Pi) - ^^^{Pl - Pi? - ie 

(28) 
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Fig. 3. Triangle diagram: Integration domain and location of singularities in the 
{w, 2;)-plane. The light shaded domains indicate the contour of integration of the 
regularized part 7,3, areas marked with (— ) enter with negative sign. The dark shaded 
area denotes critical points associated with the leading singularity of the triangle 
diagram. The two- heavy-particle (x) and the two heavy-light (•) thresholds lie on 
the boundary of the original integration domain [—1, 1] x [0, 1]. Also indicated are the 
heavy-light pseudo-thresholds (o). The corner singularities ((g)) give rise to fractional 
powers of M. 

and the integral I^b is on the other hand given by: 




It is to be noted that the integral I^a has no ultraviolet divergences. Such 
divergences can only come from the boundary contributions in the integration 
by parts over w, which in this case have to vanish because the domain of 
integration extends to the entire real axis. The low energy parts do not contain 
any fractional powers of M, as their explicit evaluation reveals the absence 
of log terms. Thus, as in the case of the bubble diagram, the separation 
into low-energy and regular parts is also possible in the triangle diagram. The 



9 



general procedure is now clear and can be easily carried over to other diagrams. 
As an illustration, we show this for the box diagram. 



2.3 The box diagram 



^ • T— ► 

Pi \ I P3 



/^+i(P3-Pl)i |fc-i(P3-Pl) 



I 
I 
I 
I 

-•- 



P2 -fc+l(P2 + P4) ^4 

Fig. 4. 

The scalar integral for the box diagram in Fig. 4 is 



/ d^ 

'4 = A* 



4-./ 



i i 
The light particle propagators can be combined to yield: 



dm? Jo J {27iY [k — q{x)]'^ — rn?{x) +ie 
1 1 



with 



q{x) = ]^{P^-Pi){l-2x) and rn{x) = - {P^ - Pif x{l - x) . 

Both ^(x) and rn{x) are small throughout the range of integration of x. Thus, 

the momentum integral is seen to be analogous to the triangle diagram. The 
procedure for separating 14 into its regular and low-energy parts, 1^ = I4 + R4, 
is therefore identical to the prescription given for the triangle diagram in the 
previous sub-section. 
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3 Discussion 



It has been shown that each basic one-loop integral / in the two-heavy-particle 
sector can be decomposed into a low-energy part / whose low-energy analytic 
structure coincides with that of /, and a regular part which consists of the 
terms in I that have factors of fractional powers of the large mass M. The 
EFTDR is now implemented by discarding R, which amounts to removing 
specific local terms in the low-energy domain. The EFTDR so defined is a 
natural regularization. 

The local terms in I can be ordered according to a power counting scheme in 
the low-energy scales generically denoted by p. For a given one-loop diagram 
the leading order term is of order p"^ with z/ = c? — 2n£ — nh + n^, where rii and 
nh are respectively the numbers of light and heavy propagators, and Un is the 
number of factors of small quantities in the numerator. This power counting 
is the standard one obtained in the 1/M expansion [11]. 

The different non-analytic terms can be classified according to their infrared 
behavior. In particular, some of them are more singular than the local terms. 
For example, the low-energy power counting for the bubble diagram would 
he 1/ — 2, but Eq. (10) shows that it is more singular than this. If one takes 
s — 4M^ to be of order p^, the bubble diagram is of order p. In the case 
of the triangle, u = 0, while /3A is more singular by one unit. As it is well 
known, these more singular behaviors are associated with the existence of the 
two-heavy-particle threshold. 

As first noticed in reference [7] , the fact that the low-energy and regular parts 
have different behaviors in M (the regular part being entirely of fractional 

power of M and the low-energy part having no such terms), implies that the 
symmetries respected by dimensional regularization are also respected by each 
part separately and, therefore, by EFTDR. 

Finally, the generalization to diagrams that involve insertions of small exter- 
nal momenta can be carried out quite easily. First, the propagators along each 
heavy particle line are separately combined using two sets of Feynman param- 
eters, and a parameter w is then used to combine the two results. Next, and 
after combining the light particle propagators, a second parameter z is used 
to combine this result with that of the heavy propagators. The integration 
over the loop momentum is then performed, and subsequently the integration 
domain in {w, z) is extended as in the case of the triangle diagram, leaving 
the integration domains of the remaining Feynman parameters unchanged. A 
more extensive presentation of this and other points in this letter will be given 
elsewhere [12]. 
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